Discrimination between evolution operators 



T. Vertesi 1 -* and R. Englman 2 ' t 

1 Institute of Nuclear Research of the Hungarian Academy of Sciences 
H-4001 Debrecen, P.O. Box 51, Hungary 
2 Soreq NRC, Yavne 81800, Israel 

Under broad conditions, evolutions due to two different Hamiltonians are shown to lead at some 
moment to orthogonal states. For two spin-| systems subject to precession by different magnetic 
fields the achievement of orthogonalization is demonstrated for every scenario but a special one. 
This discrimination between evolutions is experimentally much simpler than procedures proposed 
earlier based on either sequential or parallel application of the unknown unitaries. A lower bound 
for the orthogonalization time is proposed in terms of the properties of the two Hamiltonians. 

PACS numbers: 03.67.-a, 03.65.Ta 



I. INTRODUCTION 

Quantum information is predicated on the prepara- 
tion and identification of evolving quantum states; conse- 
quently, the determination of quantum states is of prime 
importance ([1], Chapter 6). Perfect determination is 
possible only for states that are known to be orthog- 
onal [2, 3]; in the absence of such knowledge the effi- 
ciency of imperfect determination can be improved by 
making non-projective and/or collective measurements 
whose outcomes are associated with entangled states, or 
by exploiting the quantum mechanical interference be- 
tween states [4]. Still, for perfect determination of non- 
orthogonal (pure or mixed) states an infinite number of 
state-copies are needed. Recently a probabilistic measure 
of the discriminatory process was proposed [5] . 

The above picture changes when discrimination is 
made between unitary operations; for these a finite num- 
ber of runs of the unknown gates suffice [6]- [8], a fact 
which is by itself surprising, since the number of inde- 
pendent parameters entering the unitary operations is of 
the order of the square of the number of states (assuming 
the latter to be finite). What balances this circumstance, 
is that there is a choice of states on which the unitaries 
can be tested. It now seems safe to state that while in 
certain cases the entanglement property of this chosen 
state can aid in the discrimination, perfect discrimina- 
tion is also possible without entanglement [8]. 

The present work addresses unitaries in the framework 
of the class of evolution operators and, in particular, 
those with time-independent Hamiltonians H, so that 
their form is e~ lHt . Here the time t takes over in many 
senses the role of the (discrete) number N of the uni- 
tary gates featured in the above mentioned works [6]- [8], 
but in terms of practical effort the automatically time- 
evolving unitary is much preferable to repeated subjec- 
tion to gates. The aim of this work is thus the deter- 
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mination of a minimal discrimination time between two 
different evolution operators, with liberty in the choice of 
the starting state. We identify the discrimination time 
with the time (to be denoted by t±) that two differently 
evolving states first become orthogonal. An algebraic 
treatment of the opposite issue, the recurrence of states 
under pulsed and other time dependent perturbations, 
was given in [9] and [10]. 

The relationship of this topic to time-orthogonality of 
states and to the rate of sweeping the Hilbert space is 
evident, the latter issue having already been solved for 
some states in [11]- [13] and endowed with geometrical 
interpretation in [14]- [15]. In these, lower bound expres- 
sions are given in terms of the average energy E and 
the mean deviation of the energy AE of the state. The 
states under present consideration are multi-component 
pure states, for which the existence of orthogonality times 
is studied in section III. The obvious difference (noted in 
section IV) with the previous references is that we have 
two different (in general, non-commuting) Hamiltonians, 
rather than a single one. A single-qubit state example is 
fully worked out in section V. The lower bounds tw for 
t± derived in section 6 by two different methods are sup- 
ported by the solved example (of section V). The fastest 
orthogonalizing pair of Hamiltonians turns out to be a 
rather trivial couple (section VI C). A discussion in sec- 
tion VII is mainly devoted to entanglement. 



II. TWO HAMILTONIANS 

We consider two Hamiltonians H a and that give 
rise to the evolution operators e~ A ^ and e~ B C) . Here 
A(t) = iH a t and B(t) = iH^t, where t is time. Physically, 
H a can be made different from by adding different 
excitation operators to the common system Hamiltonian 
(say, Hq). One can discriminate between the two Hamil- 
tonians (or excitation mechanisms) if, for some initial 
state ip, e~ B ^ip is orthogonal to e~ A ^ip at some or- 
thogonality time (or times) t = t± [7]- [8]. This means 

{i>\e B ^e- A ^\i>) = (il>\e iHbt ±e- iHat ^\ip) = 0. (1) 
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III. DERIVATION OF ORTHOGONALITY 



Since e A ^ and e B ^ are both unitary, their product 

(2) 



e C(t) = e B(t) e -A(t) 



is also unitary and can therefore be brought to the form 
diag(e 46,1 (*), e 1 ® 2 ^ , . . . , e^ d (*)), by a unitary transforma- 
tion matrix M{t) in which all the 9's are real (d is the 
dimension of the system's Hilbert space). When A(t) and 
B(t) are different (i.e. A{t) is not of the form: B(t)+ a 
Hilbert-space scalar, which is possibly time dependent), 
not all the 0's are the same and M(t) differs from the 
unit matrix. (In the converse case, when A(t) and B(t) 
are not different, all 6(t)'s are the same and M(i) is the 
unit matrix.) Suppose now that two #'s that differ have 
indices i and j . 

We now construct the state 



ln[e B ^e~ A ^}) that contains commutators of A(t) and 
B(t) to arbitrary order. However, to the zeroth order of 
commutators, or equivalently for the linear-in-i approx- 
imation for C(£), this quantity is a " Hamiltonian" , it 
being the difference of two Hamiltonians. To this order, 
therefore, e c ^ is an evolution operator. The speed at 
which a Hamiltonian evolution operator develops a state 
into one or several orthogonal states has been studied 
previously (e.g., [11] and references therein). One result 
is that the time t± at which an initial state (starting at 
zero energy) evolves to an orthogonal state (in units for 
which h = 1) is given by 



t± > 



2E' 



(7) 



^(t)=M + (t)v(i,0;j,a)/V2 



(3) 



in which v(i,0;j, a) is a d-dimensional column vector 
whose i'th element is e ?0 = 1 and whose j'th element 
is e ta (a being real), all other entries in the column be- 
ing zero. 

We now form the bracket 



\ e B(t) e -A(t) 



14>) = {4,\M+{t)di&g{e l9l(t \e t62{t \... 
xM{t)\il>) 



i6 d { 



iOi(t) 



(4) 
(5) 



where E is the average energy. 

The use of this (approximate) measure can be em- 
ployed in various situations, for bipartite and multipar- 
tite systems (entangled and otherwise). We shall also use 
it in the following illustrative example. 



V. ILLUSTRATION: A QUBIT 

Let us write out the Hamiltonians for the evolutions a 
and b explicitly, 



Since the 0's differ (= the cases when A(t) and B(t) are 
different), it is expected that for some t = t± they will 
differ by it and the bracket will vanish. 

The reason that the state ip has only two non-zero 
components can be explained as follows. Since 



H a = r 0a I + LJ a ^2 r ia a k , 

k—x,y,z 
k—x,y,z 



(8) 



|diag(e^«M=^ 



,i6 } (t) 



Hi 



(0) 



and the numbers IV^I 2 form a discrete probability dis- 
tribution, the set of possible values of (ip\ diag(e l6,J < - t -')|'0) 
forms a convex set in the complex plane, with the num- 
bers e t0j being extremal points. Initially, this set is the 
singleton {1}. With increasing time, the size of the set 
increases. Orthogonality is first achieved when the set 
includes the point 0. This happens when hits one of 
the edges of the convex set. The edges correspond to ip 
with at most two non-zero components. 



IV. EVOLUTION-OPERATOR ANALOGY 



Writing the (unitary) exponential e c<yt ^ introduced in 



equation (2) as e lt ^' 



-CM, 



the quantity in the square 



bracket is not time independent. This is seen by recall- 
ing the Campbell-Bakcr-Hausdorff formula for C(t) = 



Then according to Sec. II, we can write A(t) = iH a t and 
B(t) = iH t, where t denotes time. The formula which 
we intend to verify is Eq. (1), namely, for any given pair of 
H a and H we choose t = t± appropriately so that Eq. (1) 
is fulfilled. Next we omit the tq part from Eq. (8) since 
this induces merely a global phase shift in the evolution 
of the state, which has no effect on the left hand side of 
Eq. (1). 

The problem is an analog of a spin- 1/2 particle process- 
ing with frequency io a for a time t = t± in a magnetic 
field pointing in the direction r a = [r xa , r ya , r za ] (where 
the corresponding Larmor frequency uj a is proportional 
to the magnetic field strength) and subsequently process- 
ing about t = — [y'a'bi Tub , T z h] for the same time interval 
t = t± with a frequency cj&. \f a \ = = 1. The ques- 
tion is, whether the final state ipf = ip(t = 2t±) becomes 
orthogonal to the initial state ip = ip{t = 0) for a suit- 
able t±. Visualizing this situation within the Poincarc- 
spherc description, we ask whether the resulting state 
tpf points to an antipodal point on the sphere, meaning 
(V'/IV') = 0- Thus the problem can be converted to a 
problem of purely geometrical origin. 

Nevertheless, let us remain at the quantum-mechanical 
level, and express the rotation matrices exp(—iH a t) and 
cxTp(iH t) by Pauli matrices, by applying the following 
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identity: 



R P (9) 



exp(-i9fa/2) = cos(0/2)l 

-i sm(9/2)(r x a x + r y a y + r z cr z ), 



(9) 



where the operator Rp{9) in the Poincare picture defines 
a rotation by 9 about the axis f. 

After some algebra we can express the product of the 
two operators Rp b (9 b ) and R? a (—9 a ) by one operator 
R? ab {6ab) = Rr b (9b)Rr a {—9 a ) [16], where the overall an- 
gle 9 ab and the new axis f ab are given by 



cos(-0 a /2)cos(0 fc /2) 

- sin(-0 a /2) sin(0i,/2)f a • f b , (10) 



sin r a b = sin(-0 o /2) cos(0 b /2)r a 

+ cos(-0 a /2)sin(0 b /2)r 6 

- sin(-0 a /2) sin(0 6 /2)r 6 x r a .(ll) 

Let us relabel the Poincare-sphere so that r a f, points to 
the north pole. In this case, the imaginary part of equa- 
tion (9) is simply —i sin(^j L )\r ab \a z and in order to sat- 
isfy the imaginary part of (i/jf\ip) = 0, the initial state 
must take the form 



rl><*[\ T) + e ia ||>] 



(12) 



Then (ip\a z \ip) = 0. 

We now turn to the second task: to nullify the contri- 
bution from the real part of equation (9) . Let us call 7 
the angle between f a and f b . Then we have in equation 
(10) f a -f b = cos 7. With equation (10) in hand we may 
write the real part of the orthogonality condition in the 
form 

= cos(-0 a /2) cos(6> fc /2) - sin(-0 a /2) sin(0 h /2) 003(7). 

(13) 

With further algebra and observing that in our case 9 a = 
2u) a t and 9 b = 2uj b t, we arrive at the form 







cos 2 — cos(oj a 



0J b )t 



2 T 

sin - cos(w a 



LO b )t. (14) 



Let us pick from the right hand side (RHS) of this equa- 
tion the term with greater amplitude. At t = the 
RHS is equal to 1. Then if cos 2 (7/2) > sin 2 ( 7 /2), at 
t = t'™ ax = \n/(oj a — Ub)\ or else (by the converse in- 
equality) at t = t™ ax = 7r/(u> a + u>b) the expression on 
the right must be negative. Thus by continuity argument 
there must be a point (t±) between t — and t = t™ al ' 
where the RHS of Eq. (14) becomes zero. However if 



cos 2 (7/2) > sin 2 (7/2) that is 7 < ir/2 and 



, . 4-max 

is not 



becomes infinite, and the condition 
guaranteed to be fulfilled. 

This result has the following implications. If cu a 7^ Wj,, 
then the Hamiltonians a and b could be perfectly discrim- 
inated. However, if uj a = to b and the angle 7 between f a 
and f b is smaller than 7r/2 then Eq. (1) cannot be fulfilled 
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FIG. 1: Orthogonality time in a qubit against relative fre- 
quency difference r (defined in equation (19) ) of the precess- 
ing magnetic fields, for perfect alignment (7 = 0). The curve 
obtained from the solution of equation (13) is indistinguish- 
able from the evolutionary orthogonality time in the RHS of 
equation (7) . The time is in inverse units of the average 
precession radian frequencies . 



for any finite time When uj a 7^ uj b but r a = r b , the 
two Hamiltonian operators are different (by the defini- 
tion given) and they can still be perfectly distinguished 
within a finite time t, although they commute. In [7, 8] 
it is the run number N, which is broadly equivalent to 
the orthogonality time in the present work, that tends 
to infinity when the two unitaries become the same. [A 
formal expansion of equation (14) for short times gives 



t 2 = 



^a+^h- 2u> a U! b f a ■ f b 



(15) 



and shows this tendency (but the short time expansion 
is no longer valid)] 



A. Evolution time approach 

The difference Hamiltonian, which as noted previously, 
can be viewed as approximately guiding the initial state 
to an orthogonal state, is 



Hb = Kr„ - u> b f b ) 



(16) 



This can be diagonalized to give the average energy E of 
the two eigen-states (with the lower state being placed 
at zero energy) given by 



E = 



^ujI + ujI- 2uj a uj b f a ■ f b . (17) 



We first show results for a case when the evolution time 
approach is expected to give orthogonality times identical 
to the unitary operator orthogonality times. This will 
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FIG. 2: Orthogonality time in a qubit against alignment angle 
7 between the precessing magnetic fields. Dashed curve: 
solution of equation (13) . Full curve: from approximate 
evolutionary orthogonality formula equation (7) . Dotted line: 
lower bound obtained from equation (23) . The time is in 
inverse units of the average radian frequency " a 4 ^" t ' . The 
ratio of the frequencies was fixed at ^f- = 3. For negative 
angles the curves are the reflections of those shown. 



be the case when the two magnetic fields are perfectly 
aligned: 



r a ■ n = cos 7 = 1 , 7 = 0. 



The fields' Larmor frequencies still differ, 



(18) 



UJb 



= r ^ 



(19) 



so that while the a and b Hamiltonians commute, the 
evolution operators are still different. The computed or- 
thogonality times (t±), obtained from solution of equa- 
tion (14) , as function of the normalized frequency dif- 
ferences r are shown in Figure 1. When we compare 
them to the approximate (lower bound) expression given 
in the RHS of the inequality (5), using the above formula 
for E, we obtain perfect agreement (without adjustment 
of a proportionality factor). 

In the next figure (Fig. 2) we show (with the dotted 
curve) the computed orthogonality times, again obtained 
from solution of equation (14) , as a function of the align- 
ment angle 7, with the frequency difference fixed (at 
r = 0.5) and compare this to the approximate (lower 
bound) expression shown in equation (7) (the full curve 
in Fig. 2). There is still a reasonable agreement between 
the two curves. However, at non-perfect alignment an- 
gles (I7I > 0) this agreement gets spoiled for frequency 
differences smaller than r = 0.5 or — < 3. 



VI. LOWER BOUND ON THE 
ORTHOGONALITY TIME 

A. Derivation from the Anandan-Aharonov 
relation [17] 

We derive a lower bound on the orthogonality time 
t = t± for general d-dimensional systems, introduced 
in Section II. For that purpose we invoke the geometric 
quantity [17] 



s = 2 J AE(t)dt, 



(20) 



where AE(t) is the energy uncertainty in state ip defined 
by 



AE = ^h^WT^WhW 



(21) 



According to [17], the quantity s is the distance along a 
given curve C in the projective Hilbert space P as mea- 
sured by the Fubini-Study metric. In Section II, equation 
(1) can be interpreted as state ip evolving for a time in- 
terval t± under Hamiltonian H a and, following that, for 
the same interval tj_ under Hamiltonian — Hb to an or- 
thogonal state ipf. Thus one is allowed to write the right 
hand side of equation (20) as two separate time integrals, 
whose sum is 2AE a t± + 2AEbt± in our particular case. 
On the other hand, the shortest distance between orthog- 
onal states can be achieved along a geodesic curve where 
s = tt, hence in a generic evolution we have the relation 
s > tt. Summing up these facts and substituting them 
into equation (20) we obtain the lower bound 



*L > 



2{AE a + AE b ) 



tib 



(22) 



on the orthogonality time needed to discriminate per- 
fectly between two evolution operators e~ A ^ and e~ B W 
defined in Section II. 

Now let us suppose that we have knowledge only about 
the eigenvalues of the Hamiltonians H a and Hb in ques- 
tion. Having the eigenvalues of a d-dimensional Hamil- 
tonian H' (where d can not be infinite but can be ar- 
bitrarily large) gives rise to the upper bound 2AE' < 
E 'max _ E ' mm = 2uj', where E' max and E w denote, 
respectively, the highest and lowest eigenvalues of the 
Hamiltonian H' . The preceding upper bound on inser- 
tion into equation (22) yields the inequality 



t± > 



2ui a + 2ujb 



(23) 



for the time tj_ to discriminate between the pair of evolu- 
tion operators e~ A ^ and e^ B ^\ Note, that this bound 
on t± is sharp in the sense that, given uj a and LUb, we can 
always construct appropriate Hamiltonians H a and Hb 
from them so that the bound (23) for a given initial state 
ijj would be saturated. This can be achieved by choosing, 
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in particular, the Hamiltonians 

H a = Lu a \E + )(E + \-u a \E_)(E_\ , 

H b = -Lj b \E+)(E + \+u b \E_)(E_\ , (24) 

and letting the initial state \ip) = ^(\E+) +e la \E_)) for 

some a £ [0, 2ir]. We obtain for the special studied case 
d = 2 in Section V, that the optimal t± is obtained for 
the alignment angle 7 = it, a choice that indeed gener- 
ates a pair of Hamiltonians of the kind corresponding to 
equation (24). 



B. Derivation through Brody et al's formula 
[14, 15] 

Actually, in order to obtain the bound (23) we do not 
need to resort to the Anandan-Aharonov relation (20), 
as we will show in the following by using formulas from 
Refs. [14] and [15], whose derivation only requires elemen- 
tary trigonometry. The problem discussed in these works 
is the following: Consider a pair of initial and final states 
ijj] and ipF in a d-dimcnsional Hilbcrt space. The task is 
to find the Hamiltonian H on this Hilbcrt space which 
takes ipi into ipF in the shortest possible time r. It has 
been found by elementary considerations (e.g., Eq. (5) in 
Ref. ([15])), that 

T = 2arccos|(V>j|^ F )| ^ 
2tu 

where 2uj = E max — E mm denotes the difference of the 
largest and the smallest eigenvalues of H . In applying the 
above formula (25) for our situation we need to involve 
an intermediate state tp m , so that \tp m ) = e~ lHat± \ip) and 
\tpf) = e lHbt± \ip m )- Rearranging equation (25) we obtain 
two equations involving our three states, 

a a = 2arccos|(V>|V'm)| = t±2w a , 
otb = 2arccos | (V'mlV'/) I = t±2ujb ■ (26) 

If we add up the two angles, a a and a b , on the left hand 
side of equation (26), since (ip\ipf} = 0> the sum can mini- 
mally take up the value tt. Then we have t_\_(2uj a + 2ujb) = 
&a + ctb > tt- That is, t± > it j(2uj a + 2ujb)] thus we are 
back to the inequality (23) derived earlier with the aid 
of the Anandan-Aharonov relation. However, from the 
present derivation it is more transparent that in order to 
achieve the smallest value for t±, it is necessary to ex- 
press ip m , using some linear combination of the initial and 
final states tp and tp / (otherwise the angle a a + on, would 
be greater than tt). Since for a generic pair of Hamilto- 
nians H a and H b with dimensions d > 2 the condition 
that ip m is in the two-dimensional subspace spanned by 
{\ip), \ipf)} is a very severe condition to meet, the orthog- 
onality time tj_ would in general be much larger than the 
value given by the lower bound (23). 



C. Finding Hamiltonians for the lower bound 

This subsection generalizes the argument leading to 
equation (24) . 

We consider two Hamiltonians H a and Hb that give 

rise to the evolution operators e~ A ^ and e~ S( ^ . As 

before, A(t) = iH a t and B(t) = iH^t. We write 
e B(t) e -A(t) = e -c(t) = e -m c (t)t_ We assume to k now 

orthogonality lower bounds tf b = tib(E a ) for H a alone 
(and similarly i/f,(£'h) for Hb) as given in previous work 
[8]-[ll], expressed in terms of some energy E at b(e.g., the 
span of the energy spectrum, or the energy uncertainty, 
as in equation (20) ). Having the Hamiltonian H a we now 
try to find another Hamiltonian Hb such that will achieve 
the lower bound of the orthogonality time. [We shall also 
trivially find the upper bound (00) orthogonality time]. 

We recall from the Campbell-Baker-Hausdorff expan- 
sion that one can express — C(i) = \n[e B ^ tS) e~ A ^] or 
—itH c (t) = \n[e lHb e~ lHat ] in terms of a hierarchy of com- 
mutators, the zero order term being H c (t) w H a — Hb = 
H c . We also assume that H a is given and is such that 
TrH a = (which should be possible to construct for 
Hamiltonians having a finite basis) and consider only 
traceless H b {TrH b = 0). We further set Tr(H^) = 
K 2 Tr(H%), where K is an arbitrary constant. There 
may be some terms in H a and Hb that are common to 
both of them and commute with both (e.g., kinetic en- 
ergy operators that commute with spin variables, where 
only the latter constitute the differing parts in H a and 
Hb.) We discount these terms. 

Then we conjecture that Tr(H^(t)) is maximal when 
Hb = KH a , where K is negative. In equation (24) , 
K = — — . (A proof for this conjecture, provided to 
us by a referee, is given with slight modifications in the 
Appendix.) 

Then H c = (1 + \K\)H a =H a - H b , Tr{H 2 c {t)) = (1 + 
\K\) 2 Tr(H 2 ), and the lower bound for c is t% = M b {E a + 
\Eb\) = tib(E c ).(We write the modulus of E b , since only 
positive energies enter in the expressions Ub (E) , whereas 
with K negative Eb could be misinterpreted as a negative 
term.) 

Also, trivially, when K = 1, H c = and the resulting 
bra-ket is always unity, never zero. This choice (not nec- 
essarily unique) achieves the infinite orthogonality time 
upper bound. 

VII. SUMMARY AND CONCLUSIONS 

This paper treated the following: Two different Hamil- 
tonians (H a and Hb) are given and it is desired to be able 
to distinguish between them. The method of this paper 
is to form unitary operators out of these Hamiltonians 
and apply them on an initial state. Should the states 
become orthogonal, this indicates the distinctness of the 
Hamiltonians. However, two issues emerge: First , there 
is no guarantee that for any chosen initial state differ- 
ently evolving states ever become orthogonal. Secondly 
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(which has some practical implications), what are the 
times (possibly minimal times) within which one can ex- 
pect orthogonality to be obtain? By constructing some 
special states out of H a and H b , as in section III and 
shown in equation (3) , one can expect to reach orthogo- 
nality. The construction is feasible for systems with arbi- 
trary Hilbert-spacc dimensionality (and for any number 
of components) , but the achievement of orthogonality is 
not guaranteed in this work. However, the illustration 
worked out in section V for a qubit supports the ex- 
istence of orthogonality, except for effectively identical 
Hamiltonians (such as when the non-scalar parts of H a 
and H b are proportional to each other). Our formal ap- 
paratus has indicated such special cases by infinitely long 
orthogonality times, as in the text after equation (14) . 

This paper, which has quested discrimination between 
two Hamiltonians, can be regarded as an intermediate 
between past works treating state-orthogonalizations and 
unitary-discriminations. We have found ways to estab- 
lish orthogonalization times t± in general pure-state sit- 
uations and calculated t± in a single qubit model. Lower 
bounds tib of t± were proposed in terms of quantities 
arising from the two Hamiltonians. 

It is not clear whether the orthogonality times can be 
shortened by starting with entangled states (or whether 
infinite ij_'s can be reduced to finite values with entangled 
states) . From the way that maximal distinguishability is 
found for single qubit states in [6] (following cq. (12) 
there) , either through maximal entanglement or by opti- 
mal state orientation (represented there by s) , one would 
expect that the t± result found for optimal state will not 
be further reduced with entangled states. 
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where || H2 denotes the square root of the trace of the 
matrix squared (the Frobenius norm). An alternative 
statement of the theorem is, in terms of Hcrmitian A 
and B with spectrum in the half-open interval (— tt,tt], 



\n{e lA e lB )\\ 2 < || A\\ 2 + \\B\\ 2 . 



(28) 



In the notation of the conjecture, Theorem 1 states that 
for U = e tH "\ V = e- lHbt and H c (t) as above, 

\\H c (t)h = II He iHat e- iHbt )h < \\H a t\\ 2 + \\H b t\\ 2 . 

(29) 

But when H b = KH a , with K = -||ff 6 || 2 /||ffa||2, 

||-Hc(*)*||a = (1 - K)\\H a t\\ 2 = \\HJ\\ 2 + \\H b t\\ 2 , (30) 

which satisfies the equality option in equation (29) , 
showing that the choice made, H b = KH a , maximizes 
Tr(H^(t)). In the proof of the theorem (which is done 
by induction) the time t is irrelevant, so we replace the 
three Hamiltonians in the text by three related Hermitian 
matrices designated as X,Y,Z(l). These are connected 
through 



JZ{s) 



e iX e isY 



(0 < s< 1) 



(31) 



with the eigenvalues of Z(s) restricted to the half-open 
interval (— tt, tt]. In addition, Z(s) is constrained to be 
continuous over < s < 1, which requires that e lX e %sY 
has no eigenvalue equal to —1 for < s < 1, i.e., no 
eigenvalue that crosses the cut in the complex plane. (It 
may be possible to remove this constraint, but the Hamil- 
tonians in the text satisfy it.) 

We first show that for infinitesimal As, ||Z(s+As)||2 < 
\\Z(s)\\ 2 + As\\Y\\ 2 . We write 



iZ(s + As) = ln[e iZ W(l + iAsY)] 



(32) 



to first order in As. Without loss of generality we can 
consider all matrices in a basis in which Z{s) is diago- 
nal, Z(s) = diag(#i, 6 2l ■ ■ ■ j 0d)i for a <i-dimcnsional Z(s), 
with — 7r < 6j < tt. For diagonal G the (Frechet) deriva- 
tive of the matrix logarithm is given by [18] 

4- ln(G + sH) = InW (G) o H, (33) 
ass^o 



Appendix: Proof of Conjecture 

The following theorem immediately establishes the as- 
sertion in section VI C of the text that with iH c (t)t = 



\vL{e iHat e- iHht ) and the ratio Tr(H b )/Tr(H^) fixed, 



Tr(Hc(t)) is maximal when H b = KH a , where K is neg- 
ative. 

Theorem 1. Let the complex logarithm be defined on 
C\R~ , i.e., with a cut along along the negative real axis. 
For all unitaries U and V 



HUV)\\ 2 < ||ln[/|| 2 + ||lnV|| 



(27) 



where o signifies the entrywise matrix product and 
In' 1 ' (G) is the matrix of divided differences defined as 



(ln [1I (G)) jfc 



(lnW(G)) w = ! G n . 



lnG n -lnG kk ? ^ fc 



Gjj — Gkk 



(35) 



Then, putting e iZ ^ = diag(e iSl , e i02 , . . . , e i9d ), for in- 
finitesimal As, 

iZ(s + As) = ln(e iZ(s) +iAse lZ ^Y) 

= iZ{s) + iAsln [11 (e lZ(s) ) e 4Z(s) Y(36) 
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Substituting for ln [11 (e lZ(s) ) as 

0*N(e«*«))„ = OV*) 
(in^^e^W))^ = &~ i6i (37) 
we get for the two-norm 

\\Z(s + As)\\l = ^l^ + A^f 

3 

To first order in As 

\\Z(s + As)f 2 = ^|^| 2 + 2A S ^|^-| (39) 

3 3 

< ||Z(*)||| + 2A*||Z(*)|| 2 ||Y|| a (40) 

< (||Z(s)|| 2 + As||Z(s)|| 2 ||r|| 2 ) 2 ,(41) 



where equation (40) follows from equation (39) by the 
Cauchy-Schwarz inequality. Thus for infinitesimal As we 
have 



||Z(s + As)|| 2 <||Z(s)|| 2 + As||r|| 2 . (42) 

The procedure can be repeated for |jZ(s)|| 2 and so forth, 
finally giving ||Z(s)|| 2 < ||Z(0)|| 2 + s\\Y\\ 2 , so that 



||^(i)!| 2 < ||z(o)|| 2 + ||y|| 2 (43) 

= ||X|| 2 +||Y|| 2 (44) 

which was to be proved. 

(The referee pointed out that the conditions on the 
Hamiltonians being traceless are not necessary.) 
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